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Here we show theoretically that the history of solid growth during "rapid" solidification must be S-shaped, 
in accord with the constructal law of design in nature. In the beginning the rate of solidification increases 
and after reaching a maximum it decreases monotonically as the volume of solid tends toward a plateau. The 
S -history is a consequence of four configurations for the flow of heat from the solidification front to the 
subcooled surroundings, in this chronological order: solid spheres centered at nucleation sites, needles that 
invade longitudinally, radial growth by conduction, and finally radial lateral conduction to interstices that 
are warming up. The solid volume (B s ) vs time (t) is an S-curve because it is a power law of type B s ~ t n where 
the exponent n first increases and then decreases in time (n = 3/2, 2, 1, . . .). The initial portion of the S curve 
is not an exponential. 

An extremely common phenomenon in nature is the S-shaped history of areas and volumes swept by 
spreading flows and collecting flows. Examples of S-curve spreading histories are the growth of popula- 
tions 1 , the spreading of technologies 2 , and the spreading of news and information 3 , Fig. 1. The histories of 
collecting flows also exhibit S curves: examples are mining and the extraction of minerals, such as the Hubbert 
peak of oil extraction 4 . 

Rapid solidification is another common phenomenon, where the S-curve history is about the growth of solid in 
a subcooled liquid or gas. This phenomenon is classical textbook material in materials science 5 " 7 , and has 
generated a voluminous body of research dedicated to explaining the S-shaped history. The classical explanation 
is based on the Johnson-Mehl-Avrami (JMA) model of solidification 8 " 13 , which consists of postulating solidi- 
fication as a swarm of uniformly distributed spheres of solid that grow from initial nucleation sites, such that 
initially the solid volume increases exponentially. Although this model leads to a formula with two empirical 
constants that can be used to correlate experimental measurements, the S-curve remains a theoretical puzzle, 
while the theoretical basis of the sphere model is being questioned 8 : the physical meaning of the two empirical 
constants is not known. 

In this paper we go back to the state of knowledge that existed before the JMA model, and show how to predict 
the S-curve of solidification purely theoretically. The starting idea is to see solidification as a "spreading flow" like 
the flows exhibited in Fig. 1, and to recognize that this spreading flow has the same natural tendency as all the 
inanimate and animate flow systems that morph freely into configurations that provide greater access to what 
flows, over time. This natural tendency is summarized as the constructal law 1415 . Here, we ask two questions: 

First, what flow is spreading during solidification? It is not the solid, because the solid and its subcooled 
surroundings are motionless. The spreading flow is the flow of heat, which emanates from the solid surface and 
flows in all directions into the surroundings. 

Second, what is the configuration of the heat flow system? We do not postulate the configuration (e.g. 
spheres in the JMA model). Instead, we rely on the physics principle 1415 that the flow system generates and 
evolves its flow architecture in order for the solidification process to be the more "rapid" at every stage in its 
history. We show that although the solidification is triggered as small spheres around nucleation sites, at 
longer times the greater solidification rate is associated with needles and dendrites, in accord with common 
observations (e.g. Fig. 2 16 ). 

Results 

Solidification begins at point-size nucleation sites around which the solid grows as tiny spheres. As shown later in 
the Discussion section, the spherical growth slows down and is replaced by needle-shaped growth, which is faster. 
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The transition from spherical to needle-shaped growth is in accord 
with the constructal law, and defines the early part of the S. Because 
most of the S-shape of the solidification curve is due to needle solidi- 
fication, we start by examining the needle configuration. 

Needle invasion. Consider the solidification of a volume of pure 
substance that is initially in a subcooled liquid state of temperature 
T f , which is lower than the solidification temperature T s . Along the 
axis of this volume grows a solid needle with a speed (V) that 
increases monotonically with the degree of subcooling (T s — T f ), 
and which has been documented extensively (e.g. Refs. 17-21). 
The instantaneous length of the needle is x, and the length of the 
volume is L (Fig. 3). 

Because of solidification, the needle grows in length and its older 
portions become thicker. In every stationary cross-section, the 
volume of solid is isothermal at T s , and is proportional to the volume 
of liquid heated by the latent heat of solidification released on the 
solidification front, 



B 8 (t)~B f (t) 



(1) 



In other words, the latent heat of solidification (h sf ) that was released 
on the periphery of the disc of solid is now the sensible heat of the 
annulus of heated liquid, 



p 5 h sf D2 ~ PfCf (T 8 - T f )D f 2 



(2) 



The scaling relation (2) is based on the assumption that D f ^> D s , 
and it yields the constant factor in the proportionality of Eq. (1), 
namely 



D f 



■p f Cf(T 8 - T f )' 



Ps h sf 



1/2 



(3) 



The objective of this analysis is to predict the evolution of the 
solidification process, which means to predict the history of the 
volume of solid (B s ~ D^x). Because of Eq. (3), the solid volume 
is proportional to the liquid volume heated by the heat of solidifica- 
tion, (Bf ~ D 2 x). Consequently, the shape of the history function 
B s (t) is the same as the shape of Bf(t), where in accord with Eq. (3), 



Bs 
B f 



p f Cf (T, - T f ) 

Ps h sf 



= b 



(4) 



where b is constant. The time during which the liquid and solid 
thicknesses grow to the length scales D f (t) and D s (t) is 



t = 



V 



(5) 



During this time, the liquid thickness grows by thermal diffusion, 

D f ~ (a f t) 1/2 (6) 

where a f is the liquid thermal diffusivity. From Eqs. (5) and (6) 
follows the history of the volume of heated liquid, 



0! f Vt 2 



(7) 



The first conclusion is that B f and B s increase in proportion with t 2 , 
along the needle "invasion" curve shown in Fig. 4. The volume 
increase is accelerated in time, but it is not exponential. This behavior 
lasts until the needle invades the entire liquid volume, x ~ L, which 
happens during the invasion time 
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grows to the left. 




Figure 2 | Dendritic solidification oflnconel 625 in the region close to the 
melt pool 16 . 



ti 



L 

V 



(8) 



when the liquid volume length scale is Bfi ~ otf L 2 /V, with the cor- 
responding solid volume B si derived from Eq. (4). 

Consolidation. Beyond the invasion time t b the solid and liquid vo- 
lumes continue to grow in proportion with each other, but they grow 
radially laterally until D f spreads over the entire subcooled liquid. 
This is the "consolidation" process. The growth of the liquid thick- 
ness is in accord with the solution to the problem of thermal diffusion 
around a line heat source of uniform and constant temperature (the 
needle, T s ) embedded in an infinite conducting medium (Tf) 22 : 



T(r,t) - T f Ei(-r 2 /4a f t) 



T f 



Ei(-?i 2 ) 
R(t) = 2^(ot f t) 



(r > R) 



1/2 



(9) 



(10) 



where R(t) is the radius of the solid cylinder, t = 0 is the start of the 
needle solidification process, and T(r, t) is the temperature in the 
liquid. The cylindrical geometry to which Eqs. (9) and (10) refer is 
visualized in the x = constant cut through the solid shown in Fig. 3. 
Ei is the exponential integral function 



Ei(w) 



dy 



(11) 



while X is a monotonic function of the liquid subcooling Stefan 
number [Ste = c f (T s — T f )/h sf ], as shown in Fig. 5: 



r exp(r)Ei(- 



+ Ste = 0 



(12) 



The temperature distribution in the subcooled liquid is plotted in 
Fig. 5 for the case Ste = 1. The abscissa shows that the radial length 
scale (r) of the heated liquid increases as 2(oc f t) 1/2 . The volume of 
liquid heated during the consolidation increases in proportion with 
the cross-sectional area of the annulus of heated liquid, namely 
n (r 2 — R 2 ), which has the scale a f t. 

The second conclusion is that the growth of B f (or B s ) during 
consolidation process is proportional to t, which is slower than the 
t 2 growth during the invasion process. This accounts for the inflexion 
of the S curve, which is the beginning of the slowdown of solidifica- 
tion, Fig. 4. 

Plateau. The liquid volume has two dimensions, the longitudinal 
length L aligned with the needle, and the lateral (radial) dimension 
Y, across the interstices. The time when the radial dimension of the 
heated liquid reaches Y is the consolidation time t c , 

Y 2 

tc 

Otf 



(13) 



which follows from writing r ~ Y, where r ~ (a f t) 1/2 . After this event, 
the heated liquid (already Y thick) becomes warmer at a progressively 
smaller rate, as its average temperature rises to match the 
solidification temperature T s . This final arc of the S- curve of 
solidification can be predicted as follows. 



B f ,B, 




Figure 4 | The history of the volume of heated fluid during the invasion 
and consolidation processes. 
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Figure 5 | The radial distribution of temperature in the subcooled liquid, 
and the X (Ste) solution to Eq. (12). 



The heat current flows from the solidification front across the Y- 
thick liquid of average temperature T^t), and has the scale k f 7tD s (T s 
— Ti)/Y, where it is assumed that Y >> D s . The heat current is equal to 
the rate at which the energy of the liquid layer increases, namely 
pfCf7rY 2 dT 1 /dt. Integrating the resulting equation and invoking the 
initial condition T s — Ti = T s — Tfatt = ti + t c , we obtain the 
decelerating rise of the liquid temperature, 



Ti 



(t- 



- tc)/tp 



T f 



where t p is the time scale of reaching the T s plateau: 

PfCfY 3 



k f D s 



(14) 



(15) 



The heat current released at the solidification front is proportional to 
the temperature gradient across the Y-layer, therefore it has the same 
exponential decay as in Eq. (14). The same behavior belongs to the 
rate of solidification dB s /dt, and consequently the history of the solid 
volume B s (t) has the flattening tendency sketched in Fig. 4. 

Discussion 

In summary of the preceding analysis, the history of needle solid 
volume is composed of three successive periods (invasion, consol- 
idation, plateau), each with its own time scale: t b t c and t p . The 
analysis began with the needle invasion process, and it was based 
on the assumption that the solid is needle shaped, i.e. it is slender, not 
spherical. Why and under what conditions this assumption is correct 
is the foundation of the S-curve of solidification, and it is demanded 
by the constructal law. 



log(B s ) 



Plateau 



IK Consolidation 




+ t.+t r 



log(t) 



Figure 6 | The complete S curve of solidification. Note the early 
competition between the sphere and needle configurations in the pursuit 
of more rapid solidification. The last three regions (invasion, 
consolidation, plateau) are log-log representations of the S curve shown in 
Fig. 4. 

Solidification begins at point-size nucleation sites, and for this 
reason the earliest shape of the solid is spherical, of diameter D s (t). 
The sphere is surrounded by an annulus of heated liquid of outer 
diameter Df(t), which increases in proportion with (oc f t) 1/2 . Equation 
(2) is replaced by 



p s h sf D s 3 ~ p fCf (T 8 - T f )D f 3 



(16) 



and leads to the conclusion that D s (t) and Df{t) increase proportion- 
ally, and 

V-^f/ sphere \ D f/ sphere 



The solid volume B s sphere (of order D 3 ) increases as 

B s , sphere ~ b (oCf t) 3/2 



(18) 



The needle is the alternative configuration, which competes with 
the sphere. During the same time interval (0 - t), the needle would 
acquire the volume [cf. Eq. (7)] 

B s ,needle ~ OfV^b (19) 



The intersection of Eqs. (18) and (19) reveals the time scale of needle 
onset, 

OCf 
V2 



tn 



(20) 



where as shown in Fig. 6, 



Bs, sphere B s needle ? if t < t n 



B s 



needle 



> B Sj S phere j if t > t n 



(21) 



(22) 



In conclusion, in the evolution of configuration the sphere is first 
and the needle next. Not the other way around. In accord with the 
constructal law, this sequence of heat flow designs facilitates more 
effectively the flow of heat toward equilibrium, i.e. this is the evolu- 
tionary design in which the solidification is more rapid. 

The succession of four heat-flow configurations (Fig. 6) is why the 
history of the solidified volume is S shaped. The S curve is a power 
law B s ~ t n , where the exponent varies over time in this sequence: 
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n = 3/2, 2, 1 and finally n < 1. The early part of the S-curve is not an 
exponential. 
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